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1 We point out that, in analogy with spin waves in antiferromagnets, all pa- 

rameters describing the real-time propagation of soft pions at temperatures 
below the QCD chiral phase transition can be expressed in terms of static 
^NJ ■ correlators. This allows, in principle, the determination of the soft pion dis- 

■ persion relation on the lattice. Using scaling and universality arguments, we 

determine the critical behavior of the parameters of pion propagation. We 
predict that when the critical temperature is approached from below, the pole 
mass of the pion drops despite the growth of the pion screening mass. This 
fact is attributed to the decrease of the pion velocity near the phase transition. 

Introduction. — Designed for the primary goal of discovering quark-gluon plasma, exper- 
iments with collisions of large nuclei, such as those pursued at the Relativistic Heavy-Ion 
Collider (RHIC) at Brookhaven National Laboratory, also open up possibilities to study 



strongly interacting matter at extremely high temperatures and densities. In particular, it 
is hoped that at least some part of the rich phase diagram of QCD can be explored exper- 
imentally. On the theoretical side, some problems, such as finding the equation of state of 
QCD, can be effectively solved by using numerical Monte Carlo techniques. However, many 
important issues related to real-time behavior and response of high-temperature strongly 
interacting matter cannot be systematically studied by such methods. This is because lat- 
tice techniques rely on the formulation of quantum field theory in imaginary time. As a 
result, the question of what happens to the hadron spectrum of QCD, which is important, 
for example, for the understanding of certain features of the dilepton spectrum observed in 
heavy-ion collisions, cannot be easily answered in a reliable fashion. 

Fortunately, many quantities characterizing the real-time behavior of finite-temperature 
systems can be related, by exact identities, to static (thermodynamic) functions. The most 
familiar case is the relation between the velocity of sound u, the pressure p, and the energy 
density e: u = (dp/de) 1 ^ 2 . A less trivial example is that of spin waves in antiferromagnets: it 
has long been known Q that at any temperature below the phase transition, at long enough 
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wavelengths there exist low-frequency spin waves which have a linear dispersion curve, whose 
slope is given exactly in terms of static quantities. 

In this Letter, we point out that, in thermal QCD, the dispersion relation of soft pions 
can be determined entirely using static quantities. Such quantities, in principle, can be 
measured on the lattice. Using this observation, we show that the pion pole mass, which 
characterizes the propagation of the collective pion modes, decreases as one approaches the 
critical temperature, despite the well-known fact that the pion screening mass increases in 
the same limit. 

Pion dispersion from static quantities. — From the point of view of symmetry properties, 
QCD at temperatures T below or just above the temperature of the chiral phase transition 
T c is similar to a Heisenberg antiferromagnet With two light quarks (u and d), QCD 

possesses an approximate chiral SU(2)y x SU(2)a ~ 0(4) symmetry, which is broken spon- 
taneously to SU(2)y ~ 0(3) by the chiral condensate. This is similar to the 0(3) — ► 0(2) 
symmetry breaking in antiferromagnets. Moreover, the order parameter of QCD, the chi- 
ral condensate (ipip), is distinct from the conserved charges (the vector and axial isospin 
charges), which makes the real-time behavior of QCD similar to that of antiferromagnets 
(but not of ferromagnets.) 

By analogy with spin waves in antiferromagnets one can show that, at any T below 
T c , the real part of the dispersion relation of soft pions is given by 



provided the quark masses are small enough. We use the following terminology: u is the 
pion velocity (although it is the velocity only when m = 0), m is the pion screening mass, 
and the energy of a pion at p — 0, m p = um, is the pion pole mass. At zero temperature, 
u — 1, and the pole mass coincides with the screening mass. At nonzero temperature, there 
is no Lorentz invariance, and u generally differs from 1 Such pion modes with modified 
dispersion relation are termed "quasipions" in Ref. 

We shall show that the parameters u and m can be determined by measuring only static 
(zero-frequency) Euclidean correlators. In particular, m can be extracted from the long- 
distance behavior of the correlation function of the operator n a = iip^r "^, 



where x = (r, x), ip is the quark field, a, b = 1, 2, 3, r a are isospin Pauli matrices, Tr r a r b = 
2S ab , and (• • •) denotes thermal averaging. The integration over the Euclidean time variable 
t is taken in the interval (0, 1/T). In Eq. ([|) (ipip) is the chiral condensate at zero quark 



The pion velocity u is equal to the ratio of the above defined pion decay constant / and 
the axial isospin susceptibility X I5 : 




(1) 





masses. Equation (H) also provides the definition of the temperature-dependent pion decay 
constant f. 
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This is a close analog of the equation c 2 = p s /Xm for the velocity of spin waves in 
antiferromagnets. The axial isospin susceptibility X I5 can be defined as the second derivative 
of the pressure with respect to the axial isospin chemical potential [see Eq. (^) below], or, 
equivalently, via the static correlator of the axial isospin charge densities, 

S ab X I5 = JdrdV(A a (x)A b (0)), AgE^ 7 5 yf (4) 

The right hand side of Eq. (^) is free of short- distance divergences in the limit of zero quark 
masses, when Aq are densities of conserved charges. 

The derivation of Eqs. (|T])-(fD at non-zero temperature requires an analysis of the hy- 
drodynamic theory similar to the one performed in Ref . [|l[] • This approach will be presented 
elsewhere 0. In this Letter, we use an intuitively simpler (but less rigorous) derivation based 
on the effective Lagrangian approach. This approach does not allow a correct treatment of 
dissipative effects, but will be sufficient for our purpose. A somewhat similar approach has 
been used, for low temperatures, in Ref. ||. 

Derivation. — Our strategy is to first write down the most general form of the effective 
Lagrangian of pions and then relate its free parameters to the correlation functions of QCD 
by matching the partition function Z = e vv ^ T and its derivatives in the effective and mi- 
croscopic theories. The quark part of the QCD Lagrangian at finite axial isospin chemical 
potential fi I5 is given by: 

£quark = V^fB ^ ~ ($ L M$R + h.C.) + fJLj 5 A% , (5) 

where M = diag(m u ,m,i) is the quark mass matrix. The chemical potential [i I5 is coupled 
to the axial isospin charge A^ defined in Eq. (0). For simplicity, we set m u = = m q . 

We assume that, in the infrared, the pion thermal width is negligible compared to its 
energy. This has been seen in explicit calculations at low T Q . The dynamics of the pions 
is described, in this case, by some effective Lagrangian C e s, which we assume to be local, 
allowing expansion in powers of momenta. This is equivalent to the assumption that the 
correlation functions have only pole singularities, as in hydrodynamics. To lowest order, the 
Lagrangian is fixed by symmetries up to three coefficients, f t , f s , and f m , 

f2 £2 f2 

C cS = ^-TrVoSVo^ - -fTrdiW^ + ^ReTrMS , (6) 

where £ is an SU(2) matrix whose phases describe the pions. Because of the lack of Lorentz 
invariance, ff and / s 2 are independent parameters. 

The chemical potential fi I5 enters lowest-order effective Lagrangian (|6]) via the covariant 
derivative Vo in a way completely fixed by symmetries. This can be seen by promoting the 
SU(2)a symmetry in (^) to a local symmetry and treating fi I5 as the time component of the 
811(2)^ vector potential ||. The covariant derivative Vq is forced to have the form 



V £ = ftE-^ J5 (T 3 E + £T3). ( ? ) 

The structure of the Lagrangian (||) is analogous to that of the effective Lagrangian at 
finite (vector) isospin chemical potential [l 1 [0]. A significant difference between the two 
cases is that the QCD vacuum breaks the SU(2)^ (axial isospin) symmetry spontaneously. 
It is important to note, however, that the SU(2)a is a symmetry of the Lagrangian (at 
m q = 0), as good as the SU(2)y. The conservation of the axial isospin current in the 
chiral limit makes the consideration of finite /i J5 entirely legitimate. 

The pion dispersion relation following from Eq. (g) is given by Eq. d) with 



n 2 = § and m 2 = ^ . (8) 
ft f- 



Matching the second derivative of the pressure V with respect to /i J5 in QCD and in the 
effective theory, we find the relation between f t and X I5 : 

Together with the first of Eqs. (|8|) and f = f s (see below), this implies Eq. ([|). The first 
derivative with respect to m q gives 

= = (10) 

Combining with the second of Eqs. (§), we derive the generalization of the famous Gell- 
Mann-Oakes-Renner (GOR) relation to finite temperature: 

/ s W = -m 9 (#). (11) 

Finally, we need to show that / = f s . We achieve this by treating M as an external field, 
which we parameterize as M(x) = m q e lo < a ( x ) Ta . Matching derivatives of In Z, we find 

K(^^°)> = m ^ ) l i(0) =A^M^(o)); 

vr a = #7 5 r>, a (x) = ReTrzr a S(x)/2 . (12) 

Note that 7i a is defined in the microscopic theory (QCD), while a is a field of the effective 
theory. 

The correlation function of <fi a (x) can be calculated by expanding the effective Lagrangian 
in (|5[) to second order in <p a . We expect the result to match the correlator of 7r a only for 
small momenta, which means that we have to limit ourselves to zero Matsubara frequency 
and small spatial momenta, e.g., smaller than the screening mass m a of the order parameter 
a = ipip. For the static correlator of 7r a , by integrating (0) over r and using (|10D , 
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1 /■ d 3 g e i<lx 8 ab x e - m \x\ 



f%J (27r) 3 q 2 + m 2 / s 2 47r 



(13) 



We see that, by measuring the large-distance (\x\ ^> m" 1 ) static correlation function of the 
operator n a = iip'-f 5 T a ifj, we can extract two parameters of the effective Lagrangian (^): the 
screening mass m, and f s which coincides with the decay constant / defined by (|2]). The 
third parameter, f 2 , coincides with the susceptibility X I5 , which can also be expressed in 
terms of the static correlation function in Eq. (|j). From Eq. (|3]), we completely determine 
the dispersion relation of soft pions. 

Critical behavior. — For the above results to be valid, pions must be the lightest modes. In 
particular, this requires m <C m a . If m q is very small, this condition is satisfied everywhere 
below T c , except for a region very close to T c . As T — > T c from below and m a — > 0, one can 
ask the question: What is the critical behavior of the parameters u and m when T remains 
sufficiently far from T c so that the hierarchy m <C m a <C T is maintained? 

Since u and m can be related to static correlation functions, one should expect their 
critical behavior to be governed by the same static critical exponents known from the theory 
of critical phenomena. We begin by considering the critical scaling of the decay constant 
f = f s - It is defined via the behavior of a static correlator flT3|) at distances larger than 
m" 1 . In the range of momenta m \q\ m a we have [see Eq. (0)] 

[ dr alV e~ iq - x (TT a (x) n\0)) = 5 ab ^^^r . (14) 
J / q 

On the other hand, at distances short compared to the correlation length, i.e., for momenta 
such that m CT \q\ < T, the correlator of the order parameter ipip has the following scaling 
behavior: 



drdVe~ iq - x (ipip(x)ip^(0)) ~ . (15) 



We also know that in this regime the correlators of a — i/jip and ir a = iipj^T 11 ^ are degenerate, 
since they are related by the SU(2)^ symmetry, which is restored at T c . Thus the correlator 



) must match with the correlator fll4|) at the scale \q\ ~ m a . This requires 



f = A m;"(^) 2 . (16) 

The coefficient A cannot be found from scaling arguments but is finite and regular at T c . 
The exponent rj in the universality class of the 0(4) sigma- model in d = 3 dimensions (to 
which two-flavor QCD at T c belongs 0) is known: 77 w 0.03 (see, e.g., jlOlfl). 

At T — > T c the scaling laws for the inverse correlation length m a and the order parameter 
(ipip) are also known from universality, 
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m a ~ r, (17) 

W>~^, (18) 

where £ = (T c — T)/T c . Thus we find 

f „ t 2 ^ = t^ 2 >, (19) 

where in the last equation the relation 

2(3 = v{d - 2 + rj) (20) 

is used (recall that all scaling exponents can be expressed in terms of two independent ones, 
e.g., 77 and u). From this point on, we set d — 3, so / ~ t v l 2 . The is the same as the 
Josephson scaling for the superfluid density in helium |jTT|. Contrary to a naive expectation, 
/ scales differently from the order parameter {ipip) ~ t 13 . The difference is numerically small, 
due to the smallness of r\. In the 0(4) universality class in d = 3, v ~ 0.73 and (3 ~ 0.38 

Next, we point out that X Ib is finite at T = T c , where it is degenerate with the vector 
isospin susceptibility. The singular behavior of X I5 is dominated by the mixing of Aq with 
operators linear or quadratic in a or ir a in the dimensionally reduced theory describing 
infrared modes \q\ -C T. Such a mixing, however, is forbidden by the 0(4) chiral symmetry, 
as well as by charge conjugation. This is consistent with the lattice result that the vector 
isospin susceptibility is finite at T c |]I2"|| . Since f 2 = X I5 , finiteness of X I5 invalidates a 
common assumption that f t — > at T c . Note that above T c there are no propagating soft 
pion modes, so the parameters of of the effective Lagrangian (such as f t ), and hence Eq. ([5|), 
lose their meaning, even as X I5 remains well defined and finite. This is not surprising if one 
recalls that, as T — > T c from below, the domain of validity of the Lagrangian (H) (|g| m a ) 
shrinks away and disappears at T c . 

Now we are ready to find the scaling of u. Using Eq. (H), the scaling of / in Eq. (|19|), 
and the fact that X I5 is finite at T c , we find 

u 2 ~ f 2 ~ t v . (21) 

This means that the velocity of pions vanishes at T c . 

The scaling of the screening mass m can be found from the GOR relation (|11]) (recall 
that m q is assumed to be small): 

ro ^_5ffl„ m ,^. ( 22) 

In the 0(4) universality class (3 < u, which implies that the static screening pion mass grows 
(at fixed m q 7^ 0) as T —>■ T c . This fact is already known from lattice simulations of QCD. 
The pole mass of the pion, m p , scales differently: 
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m J = u 2 m 2 = „ mq t P (23) 

^75 



This means the pole mass of the pion drops as T — > T c . 

For the formulas ([2"T|)-(|2"3]) to be valid it is necessary that t <C 1. However, for any 
m q ^ 0, these formulas break down when t is so small that the condition m <C m a is 
violated. Using Eqs. ( |T7| ) and we see that this happens when t ~ m q ^ 6 or smaller. In 
the regime t -C m 1 ^ 5 the "distance" from the critical point (T = T c , m g = 0) is controlled 
by m q , but not t. The m 9 scaling of all quantities can be obtained starting from 

$ip) ~ m\ /s at t = . (24) 

Comparing this to Eq. (|18D , we see that t and m l J^ 5 have the same scaling dimension. Using 
the scaling hypothesis we can easily obtain the m q scaling by replacing t with m q ^ 5 . For 
example, 

m 2 ~ mj-^-^, ~ m^ 1 /" 5 . (25) 

(m p now has the meaning of the typical frequency of the pion mode with zero momentum. 
This mode may be overdamped in this regime.) Both masses vanish as m q — > at T = T c ; 
however, for the screening mass m 2 3> m q , while for the pole mass ml <C m q . In particular, 
near the phase transition m p m. 

The decrease of the pion pole mass may have interesting consequences for heavy-ion 
collisions. It is the pole mass of a hadron, rather than its static screening mass, that affects 
the observed spectrum. Within statistical models for hadron production, the drop in the 
pion pole mass would lead to an overpopulation of pions at low momenta, provided the 
chemical freezeout temperature T c h, at which the hadron abundances are fixed, is close 
to T c . For a crude estimate of this effect we use Aqcd ~ 200 MeV as the typical QCD 
scale in Eq. (£25|), and T ch ~ 170 MeV as an estimate for the freezeout temperature JIB] , 



which is indeed very close to T c . The shift of the pole mass near T c is approximately 
Am = m p — m n w m 7r ((m g /AQCD) 1 ^ 2 ' 5 — 1) ~ —0.3m n , where 5 ~ 5, and the pion multiplicity 
at small momenta is enhanced by roughly exp(— Am/T c h) ~ 1.3. This is a noticeable effect, 
although it is smaller than the known contribution to pion overpopulation due to the feed- 



down from the decays of resonances ||13|| . This enhancement is comparable to the effect 
of the pion chemical potential ^ ~ 50 MeV induced by pion kinetics after the chemical 
freezeout |L4] . 



Another potential consequence of the fact that pion velocity decreases at T c is the pos- 
sibility of Cherenkov radiation of pions by a hard probe moving through the hot medium 
created in a heavy-ion collision. 
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